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Abstract 

The influence of tlie gravity acceleration on the regularized energy-momentum tensor of the 
quantized electromagnetic field between two plane parallel conducting plates is derived. We use 
Fermi coordinates and work to first order in the constant acceleration parameter. A perturbative 
expansion, to this order, of the Green functions involved and of the energy-momentum tensor is 
derived by means of the covariant geodesic point splitting procedure. In correspondence to the 
Green functions satisfying mixed and gauge-invariant boundary conditions, and Ward identities, 
the energy-momentum tensor is covariantly conserved and satisfies the expected relation between 
gauge-breaking and ghost parts. A more systematic derivation is therefore obtained of the theoret- 
ical prediction according to which the Casimir device in a weak gravitational field will experience 
a tiny push in the upwards direction. 
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I. INTRODUCTION 



The use of Green-function methods in curved spacetime, and the theoretical analysis of 
the Casimir effect, are two relevant branches of modern quantum field theory. In the former 
case, some progress along the years can be outlined as follows. 

(i) In Ref. the covariant geodesic point separation is applied to evaluate the vacuum ex- 
pectation value of the ene—entu. tenso. T„„ fo. a o,a.ive scala. field in an a.bit.a.y 

gravitational field. In Ref. the same author performs regularization, renormalization 
and covariant geodesic point separation for spin 0, 1/2 and 1 fields, massive or massless, on 
an arbitrary curved background; he also finds which terms in the vacuum expectation value 
of T^i, vanish within this framework. 



(ii) In Ref. [3|, the authors obtain a momentum-space representation of the Feynman propa- 
gator G{x, x') for scalar and spin-| fields propagating in arbitrary curved spacetimes. Their 
construction uses Riemann normal coordinates with origin at the point x' and is therefore 
only valid for points x lying in a normal neighbourhood of x' . The resulting momentum-space 
representation is equivalent to the Schwinger-DeWitt proper-time representation. 

(iii) In Ref. [4i], electromagnetic and scalar fields are quantized in the region near an arbitrary 
smooth boundary. The authors find that the components of {T^y^x))^^^ generically diverge 
in a nonintegrable manner as x approaches the boundary of the manifold from the interior. 
They therefore conclude that perfect conductor boundary conditions are pathological, in 
that the distribution of energy-momentum that they would entail if they were actually to 
obtain for arbitrarily high frequencies and short wavelengths would be such as to produce 
an infinite physically observable gravitational field. 



(iv) In Ref. pi], the authors derive the symmetric Hadamard representation for scalar and 
photon Green functions, and use these representations to give a simple definition for their 
associated renormalized energy-momentum tensors. 



(v) In Ref. [6[ , the full asymptotic expansion of the Feynman photon Green function at small 
values of the world function, as well as its explicit dependence on the gauge parameter, are 
obtained without adding by hand a mass term to the DeWitt-Faddeev-Popov Lagrangian. 
Coincidence limits of second covariant derivatives of the associated Hadamard function are 
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also evaluated. 

On the other hand, an important property of quantum electrodynamics is that suitable 
differences of zero-point energies of the quantized electromagnetic field can be made finite 
and produce measurable effects such as the tiny attractive force among perfectly conducting 
parallel plates known as the Casimir effect [7|. This is a remarkable quantum mechanical 
effect that makes itself manifest on a macroscopic scale. For perfect reflectors and metals the 
Casimir force can be attractive or repulsive, depending on the geometry of the cavity, whereas 
for dielectrics in the weak-refiector approximation it is always attractive, independently of 
the geometry j^. The Casimir effect can be studied within the framework of boundary 
effects in quantum field theory, combined with zeta-function regularization or Green-function 
methods, or in more physical terms, i.e. on considering van der Waals forces ^] or scattering 
problems . Casimir energies are also relevant in the attempt of building a quantum theory 
of gravity and of the universe 

For these reasons, in Ref. we evaluated the force produced by a weak gravitational 
field on a rigid Casimir cavity. Interestingly, the resulting force was found to have opposite 
direction with respect to the gravitational acceleration; moreover, we found that the current 
experimental sensitivity of small force macroscopic detectors would make it possible, at least 
in principle, to measure such a. effect Q. Mo. precisely, the ,avitatio„a. fo.cc o. the 
Casimir cavity might be measured provided one were able to use rigid cavities and find an 
efficient force modulation method 12] ■ Rigid cavities, composed by metal layers separated 
by a dielectric layer, make it possible to reach separations as small as 5-i-lO nm and allow 
to build multi-cavity structures, made by a sequence of such alternate layers. If an efficient 
modulation method could be found, it would be possible to achieve a modulated force of 
order 10~^^ N in the earth's gravitational field. The measure of such a force, already possible 
with current small-force detectors on macroscopic bodies, rnight open the way to the first 
test of the gravitational field influence on vacuum energy 121]. In Ref. 12], calculations 
were based on simple assumptions and the result can be viewed as a reasonable first order''' 
generalization of T^y from Minkowski to curved space-time. The present paper is devoted 
to a deeper understanding and to more systematic calculations of the interaction of a weak 
gravitational field with a Casimir cavity. To first order in our approximation the former 
value of the force exerted by the field on the cavity is recovered. 

We consider a plane-parallel Casimir cavity, made of ideal metallic plates, at rest in the 
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gravitational field of the earth, with its plates lying in a horizontal plane. We evaluate the 
influence of the gravity acceleration g on the Casimir cavity but neglect any variation of the 
gravity acceleration across the cavity, and therefore we do not consider the influence of tidal 
forces. The separation a between the plates is taken to be much smaller than the extension 
of the plates, so that edge effects can be neglected. We obtain a perturbative expansion of 
the energy-momentum tensor of the electromagnetic field inside the cavity, in terms o 
small parameter e = 2ga/c^, to first order in e. For this purpose, we use a Fermi |l3, 
coordinates system (t, x, z) rigidly connected to the cavity. The construction of these 
coordinates involves only invariant quantities such as the observer's proper time, geodesic 



the 
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distances from the world-line, and components of tensors with respect to a tetrad [ij]. This 

feature makes it possible to obtain a clear identification of the various terms occurring in the 

metric. In our analysis we adopt the covariant point-splitting procedure [3] to compute 

the perturbative expansion of the relevant Green functions. Gauge invariance plays a crucial 

role and we check it up to first order by means of the Ward identity. We also evaluate the 

Casimir energy and pressure, and in this way we obtain a sound derivation of the result in 
I 

Ref . 12] , according to which the Casimir device in a weak gravitational field will experience 
a tiny push in the upwards direction. Use is here made of mixed boundary conditions on 
the potential plus Dirichlet conditions on ghost fields. Concluding remarks are presented in 
Sec. VI, while relevant details are given in the Appendices. 

II. SPIN-1 AND SPACETIME FORMALISM 

The classical action functional for the Maxwell potential A^j_{x) reads 

S[M = -] [ F^^F^^'^^d^x, (2.1) 



4 . 

where the field strength F^,^ = V f^A^ — V^A^ = A^.^—A^.,^. By virtue of the gauge invariance 
of the action, the differential operator ^/^^^ ^ is singular (see Appendix A). To take care 
of this problem one should add a gauge-breaking term, which leads to a nonsingular wave 
operator t/"^' on the potential. In the Lorenz (this is not Lorentz!) gauge [l^, and with 
the Feynman choice for the gauge parameter (see Appendix A), this is — ^(V^A^)^. Last, 



a ghost term is necessary, where x ^i^nd ip are independent ghost fields 



The full 



4 



action is therefore 
with field equations 



U'^^'As'd^x' = 0, 



U{x, x')'ijj{x')d'^x' = 
having defined (here 5"^' = g°'^6{x, x')) 



u- 



al3' 



5AJAp. 



U{x,x') = 



6^S 



-g6J'{x,x'). 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



The energy-momentum tensor is obtained from functional differentiation of the full action 
(2.2) according to 

T^'^ = (2.7) 



and turns out to be the sum of three contributions Ta,Tb, Tgh resulting from Maxwell action, 
gauge-breaking term and ghost action, respectively [2:] (no mass term occurs since we do not 
add by hand any mass for photons or ghost fields, unlike Ref. Q). This leads eventually to 
the following vacuum expectation value: 



where, on defining the photon and ghost Hadamard functions, respectively, as 



(2.8) 



H^^{x,x') = {[A^{x),A^{x')]+) = H, 



^IW 1 



jointly with 



H{x,x') = {[xix),ilj{x')]+), 



—Hj3pi-ra' + Hpr'-a/3' + -f^r/3';/3a'] , 



(2.9) 
(2.10) 



one has 



(2.11) 



(2.12) 



(TD = lim 



gap ^g,PgUr ^ ^^^^^P _ ^^^^Xp^ ^ ^ ^ ^ ^ ^ ^ ^^^^^^^^^ 



Following Christensen [2], we separate points symmetrically, i.e. 

[Aa-p, Ap.r\+ = lim ^ {[Aa';p', ^/3;r]+ + [^a;p, v4/3/;^/] + } 
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(2.13) 

(2.14) 
(2.15) 



In the applications considered in Sec. IV, the coincidence limit will be taken in two separate 
steps: first with respect to time and the coordinates spanning the plates, and eventually 
with respect to the z-coordinate (see Appendix B for details). 

In the implementation of Eq. (2.15) one needs the geodesic parallel displacement bivector 
5''^, (in general, bi-tensors behave as a tensor both at x and at x') which effects parallel 
displacement of vectors along the geodesic from x' to x. In general, it is defined by the 
differential equations 

(2.16) 



u';p 



ct(x, x') being the Ruse-Synge world function [l8j, equal to half the geodesic distance between 
X and x', jointly with the coincidence limit (boundary condition) 



lim gl, ^ [g^,] = S^. 



(2.17) 



The bivector g'^^,, when acting on a vector B'^' at x', gives the vector B^, which is obtained 
by parallel transport of B'^' to x along the geodesic connecting x and x', i.e. 



= g^,,B^' 



(2.18) 



For the reasons described in the Introduction, we use Fermi coordinates. With our choice, 
the 2;-axis coincides with the vertical upwards direction, while the (x, y) coordinates span 
the plates, whose equations are 2; = and z = a, respectively. The resulting line element 
for a nonrotating system is therefore [3] 

ds^ = -c^ (l + e-) dt^ + dx^ + dy"^ + dz^ + 0(1 x P) = r].Jx^dx'' - e-c^df, (2.19) 
V a/ a 

where the perturbation parameter e = 2ga/(?, while rj^y is the flat Minkowski metric 
diag(-l, 1,1,1). 
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III. GREEN FUNCTIONS 



For any field theory, once that the invertible differential operator Uij in the functional 
integral is given (see Eq. (A4) in Appendix A), the corresponding Green functions satisfy 
the condition 

u^jG'" = (3.1) 

and are boundary values of holomorphic functions. The choice of boundary conditions will 
determine whether we deal with advanced Green functions G~^^^, for which the integration 
contour passes below the poles of the integrand on the real axis, or retarded Green functions 
G~^^, for which the contour passes instead above all poles on the real axis, or yet other types 
of Green functions. Among these, a key role is played by the Feynman Green function G^'p ^ 
obtained by choosing a contour that passes below the poles of the integrand that lie on the 
negative real axis and above the poles on the positive real axis. If one further defines the 



Green function [1 



C^'^ = l(G'+Jfc + G'-^fc)^ (3.2) 



one finds in stationary backgrounds (for which the metric is independent of the time coor- 
dinate, so that there exists a timelike Killing vector field) that the Feynman Green function 
has a real part equal to G''^, and an imaginary part equal to the Hadamard function H^^, 
i.e. 



G'p{x,x')-G'''{x,x') . (3.3) 



This relation can be retained as a definition when the background is nonstationary; in such 
a case, however, H^^{x,x') is generally no longer real [l7l |. 

In particular, the photon Green function G\yi in a curved spacetime with metric g^y 
solves the equation [g] 

v^P/(x)G'a.' = g,,5{x,x'). (3.4) 

The wave operator results from the gauge-fixed action (2.2) with Lorenz gauge-fixing 
functional $l(^) = V'^A^, and having set to 1 the gauge parameter of the general theory, 
so that (cf. Eq. (2.5)) 

P,\^) = -^^U. + R!:{x), (3.5) 
where = g^^^V qV p{x) . Since we need the action of the gauge-field operator P^{x) on 
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the photon Green functions, it is worth noticing that 



(3.6) 



The Christoffel coefficients for our metric (2.19) read 

= Ig""' + 9Sy,p - 9p,,s) {rf'^y, + rf'^^ - ri"^S%) ■ (3.8) 

Since the connection coefficients, to first order in e, are constant, we reahze that the Ricci 
curvature tensor vanishes to this order. On expanding (this is, in general, only an asymptotic 
expansion) 

Ga.'~GS,+6G£ + 0(6^), 



we get 



so that 



and finally 



Dpxv' — dpGxv' — ^^px^iiv' — dpGxu' — ^ax^ul'^ 



QajlXu' — dadpGxu' — da 



^ PX^Hu' 



/^(O) , ^^(1) 



pM ^(0) _ p/i ^(0) 
aP^Xv',n ^ aX^iJiv',(5 



(1) pM r^(o) 

e-5^5^,dadpGfl. 



We therefore get, to first order in e. 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



where 



|—| 0/^(0) - 7(0) 

1—1 '^//i/' — -'/iiv') 



1—1 '-^/ii/' — 



with □« = r^'^^dadp = + 92 + 92 + dl 



(3.13) 
(3.14) 

(3.15) 
(3.16) 
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To fix the boundary conditions we note that, on denoting by and the tangential and 
normal components of the electric and magnetic fields, respectively, a sufficient condition to 
obtain 



0, H„ 



0, 



(3.17) 



on the boundary S of the device, is to impose Dirichlet boundary conditions on 
Ao{x) , Ai{x) , A2{x) [19] at the boundary z = 0, z = a. The boundary condition on A3 
is determined by requiring that the gauge-fixing functional, here chosen to be of the Lorenz 
type, should vanish on the boundary (see Appendix A). This implies 



(3.18) 



To first order in e, these conditions imply the following equations for Green functions: 



0, /i = 0, 1,2, W, 



0, Vz/', 



0, /i = 0,l,2, 



1 ^(0) 



(3.19) 

(3.20) 
(3.21) 
(3.22) 



hence we find that the third component of the potential A^ satisfies homogeneous Neumann 
boundary conditions to zeroth order in e and inhomogeneous boundary conditions to first 
order. 

Since jj^J, is diagonal, by virtue of the homogeneity of the boundary conditions, 
turns out to be diagonal. On the contrary, J^J, has two off-diagonal contributions: Jq^* and 
JgQ"* , so that G^^l is nondiagonal. Let us write down explicitly the expressions for the various 



components of i.e. 



„ 7(1) _ '\ _ , -^w 

""^00' ~ 2^ ' ^^00' ,00 ^ 2 00',3' 

_^(0) 
^33',0' 



.(0) 



1 



CI J d^^"^ ^ ^ ) ^^22' 00 ^ G r)()f 



00 



2^11',3' 

1 «{0) 

2 - 22',3' 

1^(0) 
2^33',3' 



(3.23) 

(3.24) 
(3.25) 

(3.26) 

(3.27) 
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'^'^30' - ~^oo',o- (3.28) 

Now we are in a position to evaluate, at least formally (see below), the solutions to zeroth 
and first order, and we get 

= V..' I ^e-(-*')+^^^-(--^^V),^^^(,, (3.29) 

having defined 

, sin K{az<) sin K{a - zy) '/ ft on\ 

gD[z,z;K,) = , 0<z,z < a, (3.30) 

K sm Ka 

, cos K(az<) cos K(a — 2;>) , 

gjsiz, z ; k) = , 0<z,z < a, (3.31) 

K sm na 

where D, N stand for homogeneous Dirichlet or Neumann boundary conditions, respectively, 

z^ (z^) are the larger (smaller) between z and z', while k± has components {kx, ky), x± has 



components {x,y), k = ^/oj'^ — fc^, and 



C = / ^e-^-(*-*')+^^^-(--^^x)$^,,, (3.32) 



diO(i^k 
(2^ 

where the $ components different from zero are written in Appendix B. A scalar field satisfies 
the same equations of the 22 component of the gauge field, hence we do not write it explicitly. 
In the following we will write simply Gy^yi and G for the Green function of the gauge and 
ghost field, respectively. 

We should stress at this stage that, in general, the integrals defining the Green functions 
are divergent. They are well defined until x ^ x' , hence we will perform all our calculations 
maintaining the points separated and only in the very end shall we take the coincidence 
limit as x' — ^ X [20]. We have decided to write the divergent terms explicitly so as to bear 
them in mind and remove them only in the final calculations by hand, instead of making 
the subtraction at an earlier stage. 

Incidentally, we note that these Green functions satisfy the Ward identity (see Appendix 

A) 

G"t';M + = 0' ^t''"' + = 0, (3.33) 

to first order in e so that, to this order, gauge invariance is explicitly preserved (the check 
being simple but nontrivial). 
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IV. ENERGY-MOMENTUM TENSOR 



By virtue of the formulae of Sec. Ill we get, from the asymptotic expansion T^ui ^ 
-^/l!!' + fl-^/li'' ^ ^(^^) (here we present only the final results for brevity, while the complete 
calculation is reproduced in Appendix B), 

<^'°'"''> - 16^ 0- ^HF^) - ^)) '^'^^(-l- 

where is the Hurwitz (^-function C,h{Xi (3) = X]^o('^ + Z^)""^- On taking the limit z' 
we get 

where the divergent term as z' ^ z can be removed by subtracting the contribution of 
infinite space without bounding surfaces [3], and in our analysis we therefore discard it here- 
after. The renormalization of the energy-momentum tensor in curved spacetime is usually 
performed by subtracting the (T^;^) constructed with an Hadamard or Schwinger-DeWitt 
two-point function up to the fourth adiabatic order . In our problem, however, as we 
work to first order in e, we are neglecting tidal forces and therefore the geometry of spacetime 
in between the plates is flat. Thus, we need only subtract the contribution to the energy 
momentum tensor that is independent of a, which is the standard subtraction in the context 
of the Casimir effect in fiat spacetime. 

In the same way (see Appendix B) we get, to first order in e: 

hm (T(^)'^^') = diag(T(i)o°,T(i)i\T(^)22^y(i)33) 

+ lim diagf-y/7r2(z-z')4,0,0,o), (4.3) 



where 



llTT^z vr cot(^)csc2(i^l 



1200 a3 3600 60 ' ' ' 
yam 7rcot(^)csc^(^) 

3600 a3 1800 a4 120 a3 ' ^ ' ^ 

2^(1)22 ^ 2.(1)11^ ^4_g^ 

(1)33 ^ {iT^ {a-2z)) 

720 a4 ■ ' 

By virtue of the Ward identities (3.33), here checked up to first order in e, the gauge- 
breaking part of the energy- momentum tensor is found to be minus the ghost part, hence 
we compute the second only, and the result is written in Appendix B. 
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V. CASIMIR ENERGY AND FORCE 



To compute the Casimir energy we must project the energy-momentum tensor along 
the unit timehke vector u with covariant components = (v^— ^oo, 0,0,0) to obtain p = 
{T>"')u^Uu, so that 

,._^,4(_^,^_if2^(M£)),o,/,. (5.1) 

720a* \ 1200 a'^ 600 a* 60 a'^ / 

where we have substituted e by its expression in terms of g. Thus, the energy stored in the 
Casimir device is found to be 

E ^ j .3E^(T-).,.. ^ -^ii (l + if) , (5.2) 

where A is the area of the plates, c/'^S is the three-volume element of an observer with 
four- velocity m^, and we have reintroduced % and c. 

In the same way, the pressure on the plates is given by 



, tt'^ he f 2 qa\ , tc"^ he f 2 qa \ , , 

P(z = = 1 + ] P(z = a) = 1 - , 5.3 

^ ^ 240 a^ V 3 c2 / ' ^ ^ 240 a4 V 3 ' ' ^ ^ 



so that a net force pointing upwards along the z-axis is obtained, in full agreement with Eq. 
(8) in Ref. [l2^, with magnitude 

F=^^. (5.4) 
720 ca3 ^ ^ 

The reader may wonder whether the pressure on the outer faces of the cavity may alter this 

result. A simple way to answer this question is to imagine that our cavity is included into a 

surrounding cavity on both sides. On denoting by b the common separation between either 

plates of the original cavity and the nearest plate of the surrounding cavity, and assuming 

that b is such that a/b << 1, but still small enough so as to obtain gb/c^ << 1, we see 

from Sec. IV that the outer pressure on both plates of the original cavity includes the 

same divergent contribution which acts from within plus a finite contribution that becomes 

negligible for a/b << 1. To sum up, the divergent contributions to the pressure from the 

inside and the outside of either plate cancel each other exactly, and one is left just with the 

finite contribution from the inside, as given in Eq. (5.3). 

As a check of the result, it can be verified that the energy-momentum tensor is covariantly 

conserved to first order in e. To this order, the covariant conservation law implies the 
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following conditions 



eO : (T(°)'^"),, = 0, 

= (^ = 0,1,2), 



e : < 



that are indeed satisfied fcf. 



1 ((T(O)OO) + (r(0)33)) + (T(l)33) 3 = Q, 



17|,|20|). 



(5.5) 
(5.6) 



VI. CONCLUDING REMARKS 



The Casimir effect for scalar fields in curved spacetime 
jy various authors, in a number of different geometries 



has been previously considered 



22 



22 



30] . More precisely, Refs. 



231] focus on a massless scalar field in half of the Einstein static universe, while Ref. 24] 



studies the same field in a Friedmann background geometry with spherical boundary, and 
the associated nonintegrable divergence in the renormalized eriergy density. Massless scalar 
fields in spherical or cylindrical shells are studied also in Refs. 25|]- 29|], with local boundary 
conditions of the Dirichlet or Robin type. Moreover, the work in Ref. 30| deals with the 
more complicated problem of a massive nonminimally coupled scalar field in between two 
infinite parallel plates moving by uniform proper acceleration. Such a scalar field is taken to 
obey Robin boundary conditions on the plates, and the interaction forces between the plates 
are investigated as functions of the proper accelerations and coefficients in the boundary 
conditions. Interestingly, for some values of these parameters the interaction forces are 
found to be repulsive at small distances and attractive at large distances 3o| . 

To the best of our knowledge, the analysis presented in this paper represents the first 
study of the energy-momentum tensor for the electromagnetic field in a Casimir cavity 
placed in a weak gravitational field. The resulting calculations are considerably harder 
than in the case of scalar fields. By using Green-function techniques, we have evaluated 
the influence of the gravity acceleration on the regularized energy-momentum tensor of the 
quantized electromagnetic field between two plane-parallel ideal metallic plates, at rest in 
the gravitational field of the earth, and lying in a horizontal plane. In particular, we have 
obtained a detailed derivation of the theoretical prediction according to which a Casimir 
device in a weak gravitational field will experience a tiny push in the upwards direction 12]. 
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This result is consistent with the picture that the negative Casimir energy in a gravitational 
field will behave like a negative mass. Our original results are relevant both for quantum 
field theory in curved space-time, and for the theoretical investigation of vacuum energy 
effects (see below). 

We stress that in our computation we do not add by hand a mass term for photons, unlike 



the work in Ref. 



2j], since this regularization procedure breaks gauge invariance even prior 



to adding a gauge-fixing term, and is therefore neither fundamental nor desirable 0, Q]. 
n agreement with the findings of Deutsch and Candelas for conformally invariant fields 
^, we find that on approaching either wall, the energy density of the electromagnetic field 
diverges as the third inverse power of the distance from the wall. It is interesting to point 
out that, in the intermediate stages of the computation, quartic divergences appear in the 
contributions from the ghost and the gauge breaking terms, which however cancel each other 
exactly. The occurrence of these higher divergences in such terms is also consistent with the 
results of Deutsch and Candelas, in view of the obvious fact that ghost fields are not ruled 
by conformally invariant operators. 



Our results, jointly with the work in Refs. [12, IMJ, are part of a research program 
aimed at studying the Casimir effect in a weak gravitational field, with possible corrections 
(albeit small) to the attractive force on the plates resulting from spacetime curvature 32] 
(cf. the recent theoretical analysis of quantum vacuum engineering propulsion in Ref. |33| ) . 
Hopefully, these efforts represent a first step towards an experimental verification of the 
validity of the Equivalence Principle for virtual photons. 
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APPENDIX A: ON WARD IDENTITIES AND BOUNDARY CONDITIONS FOR 
GAUGE FIELDS 



In the case of gauge theories there exists, on the space $ of field histories, a set of vector 
fields Qa that leave the action S invariant, i.e. 



17| 



QaS = 0. 



(Al) 



On denoting by Q*^ the components of such vector fields, and by S^i the functional derivatives 
of the action with respect to field variables v?*, Eq- (Al) implies that 



S,Q\ = 0. 



(A2) 



By virtue of Eq. (A2), the operator S^ij is not invertible, and an invertible operator Uij is 
obtained upon adding to S the gauge-breaking term 

d^x [ d*x'P''{x)ujap{x,x')P^{x'), 



(A3) 



where is the gauge-fixing functional and u;q/3(x,x') = ojafii is a nonsingular, symmetric 
continuous matrix, possibly depending on field variables to achieve full covariance of the 
formalism 



34l |. The desired invertible gauge-field operator reads therefore 



which should be considered jointly with the ghost operator 



f/" = g^p" = p", Q 



(A4) 



(A5) 



since the full action in the functional integral for the (out | in) amplitude is given by 



1 



/9 



S + -P'^uj^fiPf' + , 



with Xa and t/'^ independent ghost fields [l7| , which obey Fermi statistics if the gauge fields 
are bosonic. 

Repeated functional differentiation of Eq. (A2) yields the classical Ward identities of the 
theory, which can be used to derive remarkable identities among the gauge-field and ghost 
Green functions. For example, the first functional derivative of Eq. (A2) yields 
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(A6) 



We now re-express S^ij from Eq. (A4) and bear in mind the definition (A5) to find 

U^JQi - Ui oops P', + S,i.Q\j = 0. 



(A7) 



Restriction to the dynamical subspace, where S^i vanishes, and composition with the ghost 
Green function G"'^, for which 

f/^ G"'' = -S'^^, (A8) 

yields 

U,,Q\G''^ + u;^sP', = 0- (A9) 
Now we act on Eq. (A9) with the gauge-field Green function G^^ (see Eq. (3.1)), finding 



therefore 



3^ 



(AlO) 



Q'; = P^. G^\ 

This equation holds for all type-I theories, i.e. all gauge theories for which the Lie bracket of 
the vector fields in Eq. (Al) is a linear combination of such fields with structure constants, 
i.e. 

[Qa,Q(^]=Clf,Q„ (All) 

with C^^^ = 0. For the case of Maxwell theory in curved spacetime, the ghost Green 
function has no explicit group indices, Q''^ reduces to covariant derivatives with respect to 
the Levi-Civita connection, and cj'^ P^^G^^ yields the covariant derivative of the photon 
Green function. Thus, the Ward identity (3.33) is eventually obtained. For the sake of 
completness we write down the expansion to first order in e of the first of Eqs. (3.33): 



-"oo'.o ' -",0' 



= 0, 
0, 



TP. 



-"oo',o 



-"30',3 ^-"OO'.O "t" -",0' ~ 



< ^U^^^ = (^ = 1,2), 



r(0) 



I, -"03',0 -"33',3 2-"33' "",3' ~ ^- 



(A12) 



(A13) 
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A recipe for writing gauge-invariant boundary conditions 35| in field theory is expressed 
by the equations 

0, (A14) 



7l) (f^ 



dM 



0, 



dM 
dM = 0, 



(A15) 
(A16) 



where vr ^ is a tangential projection operator. In the case of Maxwell theory, considered in 
our paper, Eqs. (A14)-(A16) become 



[A 



kWM 



0, 



[^{A)Ui = 0, 
[A]aA/ = 0. 

On performing the familiar gauge transformation 



(A17) 

(A18) 
(A19) 



(A20) 



Eq. (A17) is preserved under (A20) if and only if A vanishes on the boundary (see (AlOV), 
since tangential derivatives and restriction to the boundary are commuting operations 35| . 
Moreover, under (A20), the gauge-fixing functional changes according to 



^{A) - = U \ 



(A21) 



where f/ is a linear differential operator. Now if A is expanded according to a complete 
orthonormal set of eigenfunctions of ?7, i.e. 

U Ui = fii Ui, (A22) 

X = J2CiUi, (A23) 

i 

the gauge invariance of the boundary condition (A18) is again guaranteed by (A19), because 

- <^CA) = U X = Y,Cifi^u,. (A24) 



The vanishing of Ui on the boundary implies therefore that both A and vanish therein, 

if was already satisfying Eg. (A18). In the full quantum theory, A should be replaced 



by two fermionic ghost fields 



13, 



19|. 
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APPENDIX B: GREEN FUNCTIONS AND ENERGY-MOMENTUM TENSORS 



1. The Fourier transform of G^^j, 



In the following, for the sake of simplicity, we always assume z' > z, hence we have to be 
careful when computing integrals and limits, but the resulting formulae become relatively 
less cumbersome. With the notation in Eq. (3.32) we find (hereafter for the sake of brevity 
we define ^ = an, s(^) = zk — ^ and s'(^) = z'k — ^) 



00' 



Qsin ^(0 



af'ijj^ cos(2^ — s — s')s^ + of'ijj^ cos(s) cos(s')s^ 



+ sin(s) a'a;' ((s' - s'') sin(2e - s') + (2^' - s'^) sin(s')) 



2 {e + a'^') {s + s') sin(0 sin(e - s') 



$03' ^ 



la 



a; sin ^{^) sin(s)((s — s') cos(^ — s') — sin(^ — s')). 



11' 



Q sin (0 

8^' 



a^oo'^ cos(2^ - s - s')s'^ - d^J^ cos(s) cos(s')s^ 



+ sin(s) l^a^a;^ ((s'^ - s^) sin(2e - s') + (s'^ - 2^') sin(s')) 
- 2{e - aW){s + s') sin(e) sin(C - s') 



$22' ^ $11' 



k30' 



2^ 



s' — s) cos(s) + sin(s) ) sin(^ — s'). 



33' _ ^(0 



(Bl) 
(B2) 



(B3) 

(B4) 

(B5) 



2 (e' - a'uj^ {s" - 1)) cos(C - s') sin(s) 
+ cos(s) ^csc(0 ((C^ - a^c^^ (s'^ - l)) cos(2^ - s') 

+ cos(s') (- (2a2a;2 + l) + a'u;' (s'^ - l) + (^^ - a^a;^) {s + s') sin(20)) 



+ 2 - a^u;^) (s + s') sin(0 sin(s') 



(B6) 
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2. The energy- momentum tensors 



In our analysis we deal with 



J V'T^) 

where the 2i factor results from the relation (3.3) between the Hadamard and Feynman 
Green functions, and T'^^^i^^' is a symmetric tensor whose components are 

,2 



= cos(e + csc(0, 

T^'^''' = -^cos(e + .-.')csc(0, 
^^(0)03' = -1 cc; csc(0 sin(e + s-s'), 

~(o),,, = _^cos(e + .-.')csc(0, 

f(o)i3' ^ _l^^sin(e + s-s') csc(0, 

= -^cos(e + .-/) csc(0, 
f(o)23' ^ sin(e + .-.0csc(O, 

= -^cos(e + «-.')csc(0. 

Now we first take the limit as t' ^ t and x'^ x±. At that stage, on taking z ^ z' 
pass the limit under the integral, then sending w — > io;, {k ^ in), and going to spherical 
coordinates: o;^ — > — cos^ ^, ky — > sin^ ^ sin^ 0, angular integration yields eventually 

(^T('^)^^)[z, z'] = -^1 fd^ex cosh(e - s - s') csch(e)diag(-l, 1, 1, -3). (B8) 
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After integrating over the ^ variables we obtain Eq. (14. ip . Analogously, upon integrating 
over the solid angle we get, to first order in momentum space. 



TT 



/.ee^smh-^(0{-3(2 + .^)e^cosh[e(.-.')]-8cosh[e(-2 + . + .)] 



60a3 

+ 8cosh[^(s + s')] 

+ 3s^cosh[^(2 + s - s')] + 32(^sinh[^(2 + s - s')] - sinh[^(s - s')] 
+ 3^s' \ - 4 sinh[^(s - s')] + 4 sinh[^(2 + s- s')] 

+ ^s'(^cosh[^(s - s')] - cosh[^(2 + s - s' 
1 



(B9) 



120a%2 



rfee'sinh-i(0 



4 (s + s') cosh(^ + s-s')+ (^~2 cosh(2 ^ - s - s') 



+ cosh(s - s') + 2 cosh(s + s')^ sinh"^(0 
+ (^~^^ + -s'^ j sinh(^ + s — s') 



(BIO) 
(Bll) 



24a%2 



rfee'sinh-i(0 



cosh(s-s') sinh-^(O) 



+ (s + s') ^4 cosh(^ + s - s') + {s - s') sinh(^ + s - s')^ 

where we have written explicitly only nonvanishing components. Integration over ^ finally 
yields Eqs. (4.4)-(4.7). 

For completeness we write also the ghost energy-momentum tensor (T^'^)gh for Dirichlet 
boundary conditions, i.e. 



(B12) 



^2^(0)00^ 
^2^(0)11^ 



gh 



(2 + cos(^)) csc^(^) 
+ ^-^^ lini 



gh 



720 24 a* 

TT^ (2 + cos(^)) csc^(^ 



720 



24 



'->z+ vr^ (z - z'Y 



— lim 



z'->z+ TT'^ (z — z' 



+ lim 



240a4 z'~>z+ TT^iz - z')^'' 



(B13) 

(B15) 
(B16) 
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and 

TT^ (3 a- 11^) TT- 



= )6QQa^;^ +40^ (a (a - 7.) - ..(a^- .) cot(-)) .«c'(-) 

+ ^ csc^(— ) C(5 + 27r'^z{a - z)) cot(— ) - 4a7r(a - 72;)] 

- lim — (B17) 
z'->z+ t:^{z- z'Y ^ ' 

/^(i)ii\ (7r2(a-2z)) 



TT 



3600 

2 



(^a(8a^ - lhaz - %z^) - '&T:z{a^ - az - z'^) cot( — ) j csc^ 



160 a° V a / a 

(lOa^ + 0^(3 + 87r^)z - 87rV(a + ^)) cot(— ) 



480 a6 

+ 2a7r(-8a^ + 15a^ + 9z^)jcsc^(^), (B18) 

{T^'^'\^ = (T«^^),,, (B19) 



647r2 



-i-j (^4a(V'(^)(l/2-^/a) - V'(^)(-l/2 + ^/a)) 
+ (a - 2^) (^V^^^) ^1/2 - ^/a) + V'^') (-1/2 + z/a)^ ^ , (B20) 
where ip^"''^{z) is the n-th derivative of the logarithmic derivative ip{z) of the F-function. 
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